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00 ' We prove that a strongly disordered two-dimensional system localizes with a localization length 

("^ ' given analytically. We get a scaling law with a critical exponent is i' = 1 in agreement with the 

("^ ' Chayes criterion v >1. The case we are considering is for off-diagonal disorder. The method we use 

04 . is a perturbation approach holding in the limit of an infinitely large perturbation as recently devised 

■4.^ ' and the Anderson model is considered with a Gaussian distribution of disorder. The localization 

O ' length diverges when energy goes to zero with a scaling law in agreement to numerical and theoretical 

^^ , expectations. 

O ■ PACS numbers: 72.15.Rn, 72.20.Ht, 71.30.+h 

pH ■ Since the discovery of the existence of a metal-insulator transition by Anderson^ originating by a significant pres- 

(— I ' ence of disorder in a sample of material, a large part of activity in condensed matter physics pertained the deep 

I . understanding of this effect originating from the localization of the conduction electrons^. 

.^H ' Several breakthroughs as the scaling theory^ appeared that improved the understanding of the metal-insulator 

■ O ] transition. Besides, weak localization theory made clear the effect of higher order corrections to the scaling theory 

-■— ' ■ itself^. More recently, the interplay between disorder and electron-electron interaction has been clarified^. 

As in all phase transitions, also this case depends critically on the dimensionality of the system. Indeed, it is 
well-known that for d < 2 the system always localizes and so it is always in an insulating phase. This result has 
been challenged by recent experimental findings^'^ showing a 2D metal-insulator transition. This effect cannot be 
fH , explained by the localization theory of Anderson and the scaling theory. A satisfactory explanation has been recently 
Q ' found numerically® through the effect of electron-electron interaction becoming increasingly important as the density 
O \ of the two-dimensional electron gas is lowered producing a rich phase diagram in this limit. 

Notwithstanding the increasing relevance that electron-electron interaction has assumed in recent years, there are a 
r/-\ \ lot of situations for disordered crystals where the Anderson model proves to be a satisfactory description of the physics 
^ . while leaving opened several problems. Current studies point toward a proper understanding of the scaling of the 
lO ' localization length with some specific parameters of the crystal as could be the energy. This scaling is characterized 
0^ , by a critical exponent v that has been obtained both by experiments^iiiS'ii and numerically^'^^'^'^'^^ giving results for 
~2 ' V ranging from 1 to 1.6 for 3D systems. 

2k \ As it is clear from the above discussion, the lack of an analytical approach prevented an analysis of the Anderson 

vQ ■ model in the regime of strong disorder. This is a rather general situation in physics. Our aim in this paper is to apply 
f~^ ' a recent devised perturbation method conceived to treat this kind of problems. This method has been initially applied 
*>••; . to quantum mechanics ^^'^^i^^ and recently to quantum field theory^*. Applications in mesoscopic physics were also 
d ' devisedi^i^Si^. 

S . The method relies on the duality principle in perturbation theory™ that, using the arbitrariness in the choice of a 

I ' perturbation, yields a dual series to the weak perturbation one. The dual series should be intended as the one holding 
^ for the development parameter going to infinity. As we showed recently, for a time-independent quantum mechanical 
problem, the dual perturbation series is given by the Wigner-Kirkwood expansionSS. We will use this latter result to 
obtain the localization length for a strongly disordered two-dimensional system in a closed form. 

The result we obtain is that the localization length scales as 1/E being E the energy and then we have z^ = 1 at 
the leading order for off-diagonal disorder. Both the localization length and the propagator are obtained analytically 
assuming Gaussian the distribution of the disorder. It is interesting to note the relevance that the Wigner-Kirkwood 
Vh ' expansion can attain beyond the already known applications in statistical mechanics. We assume this result to hold 
5^ , in two dimension. Recent numerical and theoretical results^i^'2^ point out a scaling law with a divergence at E = 
but with v ranging from 0.3 to 0.6. We note anyhow that our leading order is an Anderson model with pure disorder, 
due to the fact that we are taking a very strongly disordered system, and this difference should be expected. To 
exploit this point we compute the first order correction. In any case it should be emphasized that our result agrees 
with the Chayes criterion as it should be i/ > 1—. 

In order to make this paper self-contained we work out the duality principle in perturbation theory directly on the 
model of a single electron in a disordered material. We consider the Hamiltonian for a two dimensional electron gas 
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-— V^t(x)v*(x) + *^(x)F(x)*(x) 
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being ^(x) the field of the electron, while we take for the random potential V{x) a Gaussian distributed potential 
having 



(y(x,)) - 



and 



having the potential defined on different lattice sites. We give explicitly the probability distribution defined as 



Z(V) = Ne 
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that is defined on the sites. A^ is a normalization constant. This means that a functional integral for computing 
averages turns into an ordinary multiple integral due to the lattice. This model can be easily reduced to the Anderson 
model with off-diagonal disorder giving the same physics. 

We can give explicitly uq by the properties of the electron gas in the material by the relation^ 

(5) 
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being r the mean free time and n{E) the density of states that in our two dimensional case is n{E) = mA/nh being 
A the area and m the electron mass. 

Weak perturbation theory is straightforwardly applied to this Hamiltonian^i giving weak localization^. This can 
be obtained by writing down the Schrodingcr equation for the above Hamiltonian as 



at 2m 
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and applying the Dyson series assuming the random potential as a weak perturbation giving straightforwardly for the 
unitary evolution operator 
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being Uq the solution of the equation 
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given by 
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We recognize here the so called interaction picture. On the basis of the duality principle in perturbation theory as 
given in^^, one could ask what series is obtained when the kinetic term is formally interchanged with the random 
potential in the Dyson series. Indeed, this will give the series 
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being now Kq the solution of the equation 



in^ = l/(x)A'o 



(11) 



being 

ifo(t,io) = e-^^(")(*"*°\ (12) 

This series is dual to the Dyson seriesii and as shown in Ref;^^ it is a seniiclassical series, the Wigner-Kirkwood 
expansion, that holds as a strong couphng expansion. This series, as we showed in Refi^, has the same eigenvalue 
expansion as the WKB series as it should be. We just showed that it is dual to the Dyson series for this case. 

Our aim will be to evaluate the leading and the first orders of this series in order to have the analytical expression 
for the localization length and its first correction. Besides the physical result, our aim is also to point out another 
approach to cope with problems in the strong coupling limit that are quite common in condensed matter physics. So, 
let us give explicitly the Wigner-Kirkwood expansion till first order, one has 
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Firstly, let us consider the leading order. We take the average on the random potential to obtain 

(ifo(^,0)) = (e-^^('')*) (14) 

that has the well-known result 



{Ko{t,0))=e ^ (15) 



being 



T^ = y/AnhTTiiE) (16) 



giving a Gaussian decay in time with a time scale defined by t^ that can be seen as a coherence time. So, we can 
Fourier transform the unitary evolution operator to give 



H-oc 



{KoiE))= e'-'{Ko{t,0))dt^V^T^e-^^ (17) 



showing that the energy is strongly localized on a scale fi/r^. So, we can get immediately the localization length by 
the relation, in complete analogy with the ID case^, 

^ =- hm i-ln|(ifo(i?))P (18) 
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giving finally 

X[E)^(^Y^. (19) 
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This is the most important result of the paper showing a scaling law of the localization length with the energy 
similarly to numerical and theoretical expectation a^^i^'^'^^ and with a critical exponent f — 1. It is important to note 
the confirmation of the singularity for E = while the critical exponent in the scaling law may show differences that 
we are going to exploit computing the next order correction as at the leading order we have an Anderson model with 
pure disorder. We note anyhow that agreement with the Chayes criterion i/ > 1 is essential in our case of a very large 
disorder. We emphasize that these results apply for off-diagonal disorder. In this case the singular behavior of the 
state at E — should be expected. 

The next step is to compute the first order correction to this result through our Wigner-Kirkwood expansion. In 
order to obtain this we use the fact that the random potential is defined on a lattice and so the derivatives must be 
properly discretized. Then, we can compute the averages. So, we consider the series (fT5|) and discretize taking for the 
derivative 

^y ^ ^ V{xk,yk)-V{xk-i,yk) V{xk,yk) -V{xk,yk-i) \ ^20) 



being a the inter-particle spacing given by a = Xj ^pirn with n the density of the two dimensional electron gas. By 
noting that the only terms surviving in the averages are those in the same site, one gets the following 
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The momentum operator can be managed by a proper choice of the initial electron field granting the substitution of 
the operator itself with the eigenvalue for the successive computations. In any case we are considering the kinetic 
term everywhere negligible even if its taking into account may displace the critical point from E = after Fourier 
transform to p^ /2'm. Finally, the averaged propagator can be written down as 
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that we can rewrite, introducing the Fermi time Tp — ma^ /h, as 
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We recognize two relevant aspects of this first order correction. One has a polynomial dependence on time and then, 
at this stage, the scaling law we have obtained gets no corrections. Besides, we recognize a development parameter 
given by T^JTp that gives a working condition on the mean free time r for our expansion to hold. It critically depends 
on the inter-particle spacing. 

In order to verify our conclusion, we take the series (|23|) and put the momentum to zero (pure disorder). Then, we 
take the Fourier transform obtaining 



i,K[E)) = ^T^e 
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where we have put 



Pi{x)^--x 



(24) 
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The polynomial character of this correction changes to a logarithmic correction to the localization length that in the 
limit A — > 00 goes to zero preserving our initial result. It is interesting to point out that what we have obtained is 
also a perturbation series in the development parameter Ep/uo, being Ep the Fermi energy, and this is an asymptotic 
series in the strength of the disorder going to infinity as promised. 

In conclusion we have applied the duality principle in perturbation theory to a two dimensional gas in presence 
of a strong disorder. This has shown an application of the semiclassical Wigner-Kirkwood expansion to a strongly 
perturbed system. We have obtained a series in the development parameter Ep/uQ in agreement with expectations 
for a strongly disordered system. The localization length has been given explicitly with a power law scaling as 1/ E. 
This gives a critical exponent ly — 1 in agreement with the Chayes criterion v > 1. The method we applied is rather 
general and could be extensively applied to other problems in condensed matter physics. 
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